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Abstract 

The thermodynamic properties of Davydov-Scott monomer contacting with 
thermal bath is investigated using Lindblad open quantum system formahsm. 
The Lindblad equation is investigated through path integral method. It is 
found that the environmental effects contribute destructively to the specific 
heat, and large interaction between amide-I and amide-site is not preferred 
for a stable Davydov-Scott monomer. 
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1. Introduction 

Does the Davydov-Scott's soliton exist at biological temperature ? The 
question has attracted more interest in the last decades [J Q, S 0| • Ear- 
lier studies using finite temperature molecular dynamics showed that the 
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Davydov soliton lifetime is only few picoseconds which is too short at the 
biological temperature. The reason is the random thermal prevents Davydov 
self-trapping from occurring as, for example as discussed in jsf which showed 
that the two-quantum state might be more stable than the one-quantum 
state. Furthermore, using the standard Davydov model, some numerical cal- 
culations also indicated that soliton is stable at 310K. On the other hand, the 
analytic calculation based on trial function or perturbation methods obtained 
that soliton is stable at 300K 

The above-mentioned calculations were performed using the equilibrium 
quantum system at finite temperature [8|. The finite temperature means 
that the quantum system is in contact with environment such as thermal 
bath. The interaction of a system with its environment is given by the dis- 
sipation effect in quantum system. However, the dissipation effect leads to a 
serious problem for quantization procedure. The most appropriate theory to 
resolve this problem is the Quantum State Diffusion (QSD) based on Lind- 
blad formulation j^. The first application of QSD to the protein model has 
been done by Cuevas et.a/.j3]. Their calculation on Davydov-Scott monomer 
showed that at room temperature the semi classical approach might be a 
good approximation compared to the corresponding full quantum system. 
However the study was focused on the dynamical aspect of the system, i.e. 
the solution of Heisenberg's picture and its wave function based on the QSD 
equation. Recent studies of the anharmonic effect for the monomer has also 



been done by us using path integral and the thermodynamic function [10 
The advantage of calculating thermodynamic function is also relevant in an- 
other approaches such as the models describing the phenomena in term of 



elementary matter interactions using lagrangian It was argued that 

the thermodynamic properties should be easier to observe than another ones 
based on the wave function. Therefore, it is important to study such system 
from statistical mechanics point of view. 

The paper studies the thermodynamic properties of Davydov-Scott monomer 
including the thermal bath. The effect of thermal bath is investigated using 
the Lindblad open quantum system formalism through path integral method. 
The temperature dependencies of the system are obtained for the normal- 
ized specific heat. It is shown that the specific heat is sensitive to the size of 
coupling constant between amide-I and amide-site. 

The paper is organized as follows. In Sec. [2] the Davydov Hamiltonian 
for a one-dimensional molecular monomer is described and the thermal bath 
effect is studied in Lindblad equation. In Sec. Elthe thermodynamic prop- 
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erties are investigated using path integral method. The paper is ended by a 
short summary. 



2. Lindblad open quantum system formulation 

We use Davydov-Scott's model of the alpha-helix protein. The Davydov- 
Scott monomer is a coupled of the amide-I oscillator that expressed by the 
coordinate {x) and momentum (p) operators, and the amide-site is expressed 
by the displacement and momentum operators, Q and P, respectively. Hamil- 
tonian in the model has the form (lo| . 



1 p2 1 

H = — + -mu^x^ H h -kQ^ + yxQ , (1) 

2m 2 2M 2 ^ ^ ' ^ ^ 

where u is the intrinsic frequency of amide-1 oscillation, x is the coupling 
constant between two oscillators, m (M) is the amide-I (amide-site) mass 
and 6 is the anharmonic coefficient. Throughout the paper we also use the 
notations x' = xV'2Mf2/^ and = a/ k/M. The Hamiltonian describes the 
Davydov-Scott's monomer as a coupled harmonic oscillator. 

If the environmental or dissipation effect can not be ignored, the physical 
system is not reversible. In another words, the irreversibility implies the 
dissipation effect in a quantum system under consideration. One of the basic 
tools to introduce dissipation in quantum mechanics is the dynamical semi 
groups, and called as the Lindblad open quantum system formalism. In 
particular, the quantum system whose the wave function equation can be 
obtained from the Lindblad equation is called as QSD [9^. 

According to Lindblad formalism, the usual von Neumann-Liouville equa- 
tion is replaced by Lindblad equation or master equation in the form of, 

I = + T.(L^PL] - \l]L,P - \PL]L,) . (2) 
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p is the density function and Lj denotes the Lindblad operator which may 
neither Hermitian nor unique. In this formalism, the operator H describes 
internal dynamics, while L represents the environmental effects in the system. 

Throughout the paper, let us assume that the environmental strength of 
the amide-site is stronger than the amide-I excitation. Since L must be the 
ffist order in Q and P, we choose the Lindblad operators as follow. 
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7 is a damping parameter and v = (e 



hn/kgT 



l) ^ is Bose-Einstein distri- 



bution function with the Boltzman coefficient ks- 

The master equation in Eq. ([2]) is calculated using Feynman path integral. 
Assuming the diffusion term is dominant over the frictional damping rate 
([Q, [Q,p]] ^ [Q, [PiP\])i Eq. (|2]) is rewritten in a differential representation, 



dp 
dt 



ih 
2m 

+ 
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^2 
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imuo 
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X 
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2h 



dQ 



{xQ - x'Q') p . 
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P- 



] {x'-x'^) 



P 
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2h 2h 



{Q' - Q'')P 



(5) 



Here, the Lindblad's coefficients are 5i = j(l+2i/) / {2MQ) , 82 = 7(1 + 2u)MQ/2, 
53 = -f{l + 2u){x'^mu)/{hnf and 5^ = VmuMilx' /{h^). The propagator of 
Eq. ([5]) is given by. 



K{x,x';Q,Q') 




V[x]V[Q]exp 



— / t I -mx^ — —muj'^x^ 



2 

■x'xQ 



(6) 



where u'^ = u'^ + iSs/imh), = + i62/{Mh), M = M + ih^/6i and 
X = X + 64/2. Making use of the Gaussian approximation, onl y th e classical 
path of amide-site (Q) contributes to the interaction term [l2|, [l3|. It yields, 

K{x,x';Q,Q') = K^Kq 



X 



P[x]exp 
I?[Q]exp 



-mx 



-U X 



(7) 



The pr opag ator Kq is just a harmonic oscillator, and the solution is well 
known 13|, [ij] , 



xexp 



iMQ 



2n\sm{nt)\ 



iQ'^ + Q^) cos{Qt) - 2QQ' 



where ||a;|| denotes the largest integer smaller than x. 

The propagator of is a driven harmonic oscillator and the solution 
is also known [l^ 14 1. The driven function Q is the classical solution of 
equation of motion (EOM) of Q, that is a harmonic oscillator. Taking the 
solution Q = Qosm{Qt) and substituting it into the path integral solution of 
a driven harmonic oscillator, then the propagator becomes 13|, Il4j |. 



K{x, t; x', 0) = exp 



TT A 


ut 
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2TTh\ sin(a;t)| 



(9) 



where Sc\ is given by, 

muj 



2 sin(wt) 

xQox' 



[(x^ + x'^) cos(ci;t) — 2xx'~\ 

ujcos{ujt) sin(nt) — flsm{ut) cos(fit) 
flsm{ujt) — ujsm{Qt) 



{Q'^ - ^2) sin(wt) 

xQox 



(10) 



x'Ql 



'mu{Q'^ — u"^) sin(a;t) 



A cos(wt) sm\nt) + B sm{ujt) sin(2fit) 1 
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with A = + l)/(4(f22 - cl;^)) and B = + u^)/{8n{n^ - u^)). 

Combining Eqs. and (HU]) . the propagator K{x,x';Q,Q';t,t') be- 
comes, 



K{Q,x,t;Q',x',0) = exp 



—tTT 



1 1 

2 + 2 



TT 
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xexp 
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3. Thermodynamic properties 



The discussion of a system interacting with heat bath is characterized by 
the temperature T. The state of those systems is therefore given by an equi- 
librium density matrix which can be obtained by performin g a t ransformation 
t — )■ r = —ih(3 in the propagator with (3 = l/{kBT) [isl, 14|. The density 
matrix is actually the propagator with p(x, x, Q, Q, r, 0) = K{Q, x; Q, x, r, 0). 
Substituting t ^ r = —ih/3 into Eq. (ITT]) , one obtains, 



{x,Q) 



exp 



1 i 

-ITT I 

2 2 



TT 



2nhsmh{nhl3) \/ 27ihsmh{uh(3) 
xexp <^ — taxi\i{-Vthl3)Q'^ — 



tanh(-a;/i/3)a;^ 



+ 



n ?,m\i^{]ph(5) - Q sinh(fi;i/3) sinh(^cu/i/3) 



a; 



+ 



x'Ql 



uh/3 
4sinh(a)/?./3) 



An smh^{nhf3) coth{uh(3) + B iwh{2nni3) 



(12) 



Having the density matrix at hand, in the statistical mechanics one can 
consider partition function 12| . 



(13) 



Substituting Eq. ( IT^ into Eq. ( IT5]) and using the hyperbolic manipulation 
yield. 



1 



1 



2sinh(i^]n/3)2sinh(|w;^/3) 



■ exp 
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+ 



mujh?{p? — Cj"^) sinh(a)/i/3) 



by borrowing the Gaussian integral, j dxe'XY>{—ax^+hx+c) = ja exp [6^ / (4a) + 
c]. 

For an open quantum system, such as the Davydov-Scott monomer, the 
changes of the surrounding environment entropy must be taken into account. 
The effect is conveniently incorporated by the specific heat which is also 
experimentally measurable. It describes the quantities of heat that must be 



added to a system in order to increase its temperature, and defined as [12 

InZ 



C = keP' 



a/32 



Bringing Eq. ( !T4|) . one gets, 

C_ P 



+ 



w 



sinh^ / sinh^(w) 



+ rm 



(15) 



(16) 



where / = flhl3/2, w = ijjhf3/2 and. 



m 



2h dl3^ 



sinh^ I — Co sinh(2/) sinh w 
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+ 



AVL sinh2(2/) cosh(2u;) + B sinh(4/) sinh(2w) - \w 
mCbh(p?' — a)^) sinh(2w) 



;i7) 



For small coupling case we have approximately sinh(x) ~ x, cosh(a;) ~ 1 and 
tanh(x) ~ a; to yield. 



M ( _ 3^2-2 _ 3~4^ I ^3 



The specific heat in Eq. ( fT6l) can be read as, 

Camide— site ~l~ C*a,mide— I ~l~ C*] 



mixing • 



(18) 



(19) 



If there is no coupling between amide-site and amide-I, i.e. x = 0, this is 
just the total specific heat of two independent harmonic oscillators. 
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Figure 1: The temperature dependence of normalized specific heat for the first (left), 
second (middle) and total (right) terms in Eq. ([T^ for various values of 7. 



The result is depicted in Fig. [T] for three cases corresponding to the 
values of 7, that is 7 = 0, 7 < and 7 > ^2. The case of non-zero but 
small 7 has similar behavior as non-damping case. All of them coincide 
each other at low temperature and tend to be asymptotically constant at 
high temperature. On the other hand, large 7 reduces the specific heat 
for the whole region of temperature. This means the environtmental effect 
contributes destructively to the specific heat, and the system requires less 
energy to increase the temperature to reach the equilibrium. 

The application of the present model to the a— helix protein requires 
proper knowledge on the coupling constant between amide-I and amide- 
site [3|]. There are some attempts to determine its allowed range through 
several methods like the Ab initio calculation and also the extraction from 
the experimental data as well. The value is found to be within 7 pN and 
62 pN Considering the amide-site mass 5.7 x 10~'^^ kg and the string 
constant k = 58.5 Nm~^, one immediately obtains the coupling constant 
x' = 1.7,11.6,18.2 corresponding to x = 0.5,2.5,5.5 pN. Note that most 
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of previous works takes x = 62 pN. However, those works do not take into 
account the thermal bath effect. On the other hand, including the thermal 
bath as done in the present paper enhance the contribution of amide-I and 
amide-site interaction. 

The result is depicted in Fig. [1] showing the temperature dependencies 
of the normalized specific heat for various values of 7. The left, middle and 
right figures correspond to the contribution of the first term, the second term 
and the total in Eq. ( [T9|) respectively. From the figures one can conclude that 
the results are sensitive to the size of coupling constant x intermediate 
temperature. Especially, the amide-I is more affected than amide-site. The 
reason is because amide-I has higher frequency than amide-site. Moreover, 
amide-I is also suppressed significantly by thermal bath contribution which 
indicates the dependencies of system frequency on the effect of thermal bath 
as already pointed out by Ingold et. a/.through Caldiora-Lenggets formalism 
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It should be remarked that at low temperature region large environmen- 
tal effect induces an anomaly, that is the specific heat is getting negative. 



This anomaly has also been observed by Ingold ei.a/.[l5| for free harmonic 



oscillator using Caldiora-Lenggets formalism, and by us using full-quantum 



approach and the Lindblad formulation of master equation [10 



4. Summary 

The interaction of Davydov-Scott monomer with thermal bath is inves- 
tigated using the Lindblad open quantum system formalism. In contrast 
with previous work by Cuevas et.a/.0], the statistical partition function is 
calculated instead of solving the EOM itself. 

Using path integral one can calculate the propagator of the Lindblad 
equation. Under an assumption that the diffusion term is dominant, we have 
shown that the environment contributions shift the kinetic, potential and 
interaction terms in the lagrangian. The mixing term is survive only if the 
frictional damping rate is taken into account. 

In the open quantum system the damping coefficient 7 represents the 
relaxation time due to interaction with the environment. Non-zero 7 con- 
tributes destructively to the specific heat. The higher value of 7 corresponds 
to the shorter relaxation time, and it induces specific heat anomaly as pointed 



out in previous works [15|, llO[. It is also found that large interaction between 



amide-I and amide-site is not preferred for a stable Davydov-Scott monomer. 
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